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Spin ice illustrates many unusual magnetic properties, including zero point 
entropy, emergent monopoles and a quasi liquid-gas transition. To reveal the 
quantum spin dynamics that underpin these phenomena is an experimental chal- 
lenge. Here we show how crucial information is contained in the frequency de- 
pendence of the magnetic susceptibility and in its high frequency or adiabatic 
limit. These measures indicate that monopole diffusion is strictly Brownian but 
is underpinned by spin tunnelling and is influenced by collective monopole in- 
teractions. We also flnd evidence of driven monopole plasma oscillations in weak 
applied field, and unconventional critical behaviour in strong applied field. Our 
results resolve contradictions in the present understanding of spin ice, reveal 
unexpected physics and establish adiabatic susceptibility as a revealing charac- 
teristic of exotic spin systems. 



1 Introduction 



In spin ice materials like Ho2Ti207 or Dy2Ti207^~^ magnetic rare earth ions (e.g. Ho, Dy) 
occupy a lattice of corner-linked tetrahedra. In the low temperature spin ice state two 
atomic magnetic moments or 'spins' point into, and two point out of each tetrahedron (Fig. 
1). This is equivalent to the ice rule that determines proton configurations in water ice^'^, 
and hence spin ice has a residual entropy equal to the Pauling entropy of water ice^. The 
thermodynamic properties of spin ice are well described by a classical spin Hamiltonian with 
a dominant dipole-dipole interaction^'^. The self-screening of the latter establishes the ice 
rule ground state but this property does not extend to excited states^. A spin flip out of 
the ice rule manifold creates a dipolar magnetic excitation that may fractionalise to produce 
free defects (Fig. 1). These inhabit the diamond lattice formed by tetrahedron centres and 
behave as magnetic monopoles on account of the integrated dipole-dipole interaction^. 

The spin ices are part of the family of rare earth pyrochlores, a series of frustrated magnets 
for which collective quantum effects have been widely discussed Recent theoretical 
16,17 ^Qgg j^Q^ possibility that such effects may be relevant to Ho2Ti207 and 

Dy2Ti207, but to a good approximation, the monopoles may be treated as classical objects, 
with local quantum mechanics setting local parameters such as attempt frequencies. 

The magnetic monopole current density in spin ice is defined as the rate of change of 
magnetization: J = 5M/5t, with the conductivity proportional to the monopole density^. 
However, even in an infinite system, magnetic monopoles in spin ice cannot sustain a direct 
current, on account of the destruction of the spin ice entropy by magnetization of the system^. 
This means that direct current (dc) 'magnetricity' in spin ice^^'^^ is necessarily transient 
Alternating current (ac) magnetricity does not suffer from this limitation as monopoles can 
in principle be driven indefinitely back and forth by an oscillating magnetic field. The theory 
of ac-current^ has not yet been tested as existing ac-magnetization studies either precede 
the theory ^^'^^ or focus on the low temperature regime ^^'^^ where complicating factors are 
expected In Section 2 we present the first experimental test of the theory of Ref.^, where 
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we rigorously confirm a number of ideas and arguments about monopole diflFusion^^'^^'^^ and 
spin tunnelling and derive new information on the microscopic processes involved. 

As shown in Fig. 2, our dynamical magnetization measurements also estimate the isother- 
mal susceptibility Xt and the adiabatic susceptibility xs- While the former is a much 
discussed magnetic response function, the latter is typically neglected. Nevertheless, our 
experimental data clearly show (Fig. 2) that xs is finite, with the ratio Xs/xt increasing 
with increasing applied static magnetic field. In Section 3 we report a striking contrast 
between the temperature dependence of xs and Xt in weak applied field, showing that Xt 
is best interpreted as a spin response, while xs is best interpreted as a monopole response. 
This contrast has its root in the fact that configurational entropy ultimately confines the 
monopoles when they are driven by a magnetic field ^. 

In strong applied field along the cubic [111] direction, spin ice exhibits a liquid-gas type 
phase transition with a critical point at /io^c = 0.929 T, Tc = 0.36 K. This transition 
has been interpreted as a monopole condensation^ and has been treated in renormalisation 
group theory^^. In Section 4 we extend our comparison of Xt and xs to the 'supercritical 
regime' at T > T^, where we observe strong signatures of critical behaviour and find that 
monopoles behave increasingly like dipole pairs, in agreement with comments of Ref.^^. 

The temperature and field regimes probed in this paper are illustrated in Fig. 1, where we 
broadly define regimes of monopolar response and dipolar response. It should be emphasised 
that, as in other cases where novel quasiparticles accurately account for the low-energy 
physics, monopole and spin descriptions are never in conflict. Instead, certain properties are 
best discussed in terms of spins and others are best discussed in terms of monopoles. One 
result of our study is to clarify how this division should be made. 

2 Magnetic Relaxation in Zero DC-Field 

According to Ref.^ (see also Ref.^^) the frequency {u) dependent susceptibility arising from 
monopole currents should be described by the following equation with a — xs — ^ and r a 
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relaxation time: 

X(^) - Xs 1 



(1) 



Xt-Xs 1 + (icjr)i-^* 
With finite xs and parameter a this expression coincides with those of Casimir - Du-Pre^^ 
and Cole - Cole^^ for ordinary magnetic relaxation, but there is an important difference: 
here it is simplest to understand the behaviour of r in terms of monopole dynamics rather 
than interacting dipoles. In particular the relaxation rate may be written — fiQuQx/Vo 
where u is the monopole mobility, Q the monopole charge, x the total monopole density 
per diamond lattice site, and Vq the volume per lattice site^^. The density x{T) evolves 
with temperature in a way that cannot be expressed in closed form^^. However, we could 
measure x{T) by fitting specific heat data to Debye-Hiickel theory^^'^^ so that by dividing our 
measured r by x(T) we were able to derive the mobility u{T) as a function of temperature 
(Methods). Representative experimental data and results, along with characteristic fits, are 
shown in Fig. 3. 

Referring to Fig. 3, the model fits the experimental data well at T > 3.5 K but describes 
only the high frequency part at the lowest temperatures (Supplementary Information, Fig. 
SI). At T > 10 K the apparent mobility diverges in accord with an expected Orbach 
type spin flip process ^^'^'''^^ that is not considered further here. At lower temperatures 
u(T) becomes accurately proportional to 1/T which is consistent with the Nernst-Einstein- 
Smoluchowski equation for the Brownian diffusion of magnetic monopoles: 

DQ ... 

Here, the diffusion constant D is temperature-independent, as shown in the inset of Fig. 
3. The athermal diffusion constant shows that the observed temperature dependence of 
the magnetic relaxation in this temperature range, is completely accounted for by the 
temperature evolution of the monopole density and the temperature factor characteristic of a 
diffusive process (Supplementary Information, Fig. S2). This general behaviour is insensitive 
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to small applied field (typically ij^qH < 50 mT, Fig. 3). Writing D = v^o? /Q where a is 
the diamond lattice constant and z/q the monopole hop rate^^'^^, we find a temperature 
independent hop rate of z/q = 2.43(6) x 10^ s~^. This athermal hop rate may be treated as 
evidence of quantum tunnelling of the spin involved in the monopole hop (Fig. 1). 

Our results are fully consistent with the theory of Ref.^ and the numerical analysis of 
Ref.^^ but also indicate an essential refinement that must be made. That is, we find a finite a, 
suggesting a significant dispersion of relaxation times, as previously observed rather than 
the single relaxation time assumed in the theory. However, the theory neglects monopole 
interactions which might be expected to infiuence the hopping rate of individual monopoles. 
To capture this we assume that in zero applied field, spins are fiipped by transverse fields 
arising from the dense ensemble of atomic dipoles, and we decompose the instantaneous local 
transverse dipolar field as follows: 



where is an effective field that causes fiipping at rate z/q^ determines a static dispersion 
in the latter field, and h^^fx determines the local field arising from the monopole gas, which 
at the level of Debye-Hiickel theory scales as ^Jx. Assuming uncorrelated contributions it 
may be shown (Methods) that the variance of the logarithmic relaxation time is given by: 



where the subscripts 1 and 2 refer to the fields h\ and /i2, respectively. Our measured Q^(T) 
may be transformed^^ to give the quantity on the left (Methods) and hence we can test 
the above expression. Fig. 4 confirms a very satisfactory agreement between theory and 
experiment in zero and weak applied field, with the fitted G2 increasing rapidly in an applied 
dc field of 1 mT, but thereafter more slowly. Possible ambiguities in our interpretation are 
discussed in the Supplementary Information (Fig. S3), but at the very least we may conclude 



H = Ho{l + hi + h2^/x) 



(3) 



(4) 
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that observed dispersion of rates is in large part a monopole property. 



3 Isothermal and Adiabatic Susceptibilities 

The isothermal susceptibility Xt extracted from the fits to theory is in close agreement with 
the directly measured Xt (Supplementary Information, Fig. S4). In theory the isothermal 
susceptibility is twice the Curie susceptibility^, Xt = 2C /T, but recent work^^ has estab- 
lished that in spin ice there is a crossover from a Curie constant C at very high temperature 
to the expected 2C at low temperature. Our results indicate that C ^ 4.25 in the tempera- 
ture range explored (Supplementary Information, Fig. S4) which we interpret as evidence for 
this crossover: however a much more detailed study of the Curie Law crossover in Dy2Ti207 
would be worthwhile. The Curie like Xt is of course characteristic of a spin system: indeed 
there is no direct monopole signature in this quantity. This may be traced to the configura- 
tional entropy in the problem, which in applied field confines the monopoles^'^^, making the 
magnetic response spin-like at long time. 

The thermodynamic adiabatic susceptibility xs is the ac-susceptibility extrapolated to 
infinite frequency (or more strictly to a frequency where spin-spin relaxation is active but 
where spin-lattice relaxation is not^^). Fig. 5 (bottom) illustrates a striking correlation 
between our measured adiabatic susceptibility Xs{T) and the measured monopole density 
x{T) (Supplementary Information, Fig. S5 and S6, for further discussion). Thus we find xs — 
Xox(T) with xo — 0.030(1), a temperature-independent constant. In a monopole picture, 
we may imagine a frictionless, and hence reversible, displacement of magnetic monopoles 
by distance r in the applied field - like a driven plasma oscillation. If we write the force 
on a positive monopole as iiqH{uj)Q — Kr{uj)^ where K is the force constant, then the 
magnetization is M{uj) = {x/Vo)Qr{uj), from which xs = x/hqQ'^ /KVq, as observed. From 
the value of xs^ we find K ^ 0.12 Nm~^^ implying an energy barrier between lattice sites 
at a distance r = a/2 of order 200 K. The latter seems too large to be a Coulombic barrier, 
and is more likely connected with the crystal field energy scale of several hundred kelvin. 
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The frictionless oscillation of the monopole ensemble is reminiscent of a plasma oscillation 
in an electrical plasma, though the absence of an accelerative term in the equation of motion 
means that the monopole plasma oscillation cannot occur in the absence of a driving field. Of 
course a finite xs ^ magnetic system can always be formally represented as an oscillation 
of magnetic charge, but in this case our result shows it to be associated with the motion of 
recognisable positive and negative magnetic monopoles. 

In a magnetic system xs is always less than the isothermal susceptibility Xr, as it obeys 
the thermodynamic relation: 



where Ch > is the specific heat at fixed applied field H. For a paramagnetic rare earth salt 
a typical behaviour of xs would be to roughly track the increase of Xt as T — )► according 
to the Curie Law Xt = C/T. The striking difference we observe between Xt{T) and Xs{T) 
(Fig. 5) refiects transition from spins to monopoles as the natural variables by which to 
describe the magnetic response, monopoles being more appropriate at high frequency. 

Nevertheless, we can explore the origin of xs in spin language if we consider a monopole 
as a label for a set of 'fiippable' spins (Fig. 1). We assume that the adiabatic susceptibility 
is equal to the isolated susceptibility, which in a semi-classical approximation is given by^^: 



Here Z is the partition function, n labels the energy states of the system, and is the 
magnetic moment per unit volume of the state n. If the ground state is assigned null 
moment and the monopole excited state is assigned dM^/dH = xo? where Xo is temperature 
independent, then we obtain our experimental observation that xs — Xo^iT). Since — 
V~^dEn/d{fioH) (where V is volume) our result reveals a quadratic term in the energy per 
monopole: E'^ = (l//ioXo/2)i?^. 



Xs^Xt- 



T{dM/dTf) 
Ch 



(5) 




(6) 
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A quadratic energy expression generally indicates 'stretchable' magnetic moments. A 
small quadratic (Van Vleck) term is expected for a free Dy^+ ion through mixing of the 
ground state with states of higher total angular momentum J. However, in our case, the 
observation that only flippable spins contribute to Xo and that flippable spins and non- 
flippable ones are distinguished only by a thermal energy scale at these temperatures, appears 
to rule out any single spin mechanism. It is interesting to note that the monopole spin texture 
is predicted to produce an electric dipole^^ and it appears from our result that it is associated 
with a magnetic polarisability as well. 

These findings have very important ramifications for the monopole description of spin 
ice. Previous work by neutron spin echo^'''^^ on Ho2Ti207 and /iSR^^ on Dy2Ti207 has 
suggested a high frequency response that thermally evolves to low temperature, and that at 
first sight seems disconnected from the monopole picture. However, our results indicate that 
the dynamical spectrum in the approach to the high frequency limit is fully accounted for by 
magnetic monopoles, and they clearly explain the thermal evolution observed in the previous 
work. A very recent thermal conductivity study indirectly estimates a diffusion constant 
for magnetic monopoles that is much faster than ours, but our adiabatic susceptibility results 
show that there is no necessary contradiction, as monopoles mediate a dynamical response 
over a very broad frequency range. Finally our results rule out any significant spectral weight 
beyond that associated with monopoles, contrary to a recent proposal 



4 Adiabatic Susceptibility in Applied Field 

Fig. 1 shows the spin ice phase diagram for a dc-magnetic field applied along the cubic [111] 
direction. A small applied field orders one spin per tetrahedron in the pyrochlore structure, 
but maintains the ice rule of two spins in and two out per tetrahedron, thus creating the so 
called 'kagome ice' phase of two-dimensional disordered sheets, which still possess residual 
entropy With increasing field at T = there is a breaking of the ice rules, pictured 
as the flipping of one spin per tetrahedron, to create an ordered 'three in, one out' state. 
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Extending from this point a first order phase transitions that terminates in a critical end 
point. The positive slope of this fine reflects the destruction of the spin ice entropy by 
the applied field, according to the Clapeyron equation. In the monopole representation, 
the applied field tunes the chemical potential of monopole-antimonopole pairs such that the 
increased monopole density drives a first order condensation from a sparse monopole fiuid 
to a dense 'liquid' (or perhaps better, ionic crystal) of alternating positive and negative 
monopoles^. The detailed theory of magnetic relaxation near the critical point predicts 
mean field critical exponents modified by logarithmic corrections. Here we are interested 
in the supercritical region at temperatures well above the critical point, where the system 
may be described as a dense monopole plasma. Recently, a peak in the ac-susceptibility at 
finite frequency was observed in this region (Matthews, M. & Schiffer P. unpubhshed). We 
examined the behaviour of the adiabatic susceptibility as a function of field in this regime, 
to compare it with our zero field measurement. 

In weak fields {/noH < 0.3 T) the thermal evolution of Xs{T) shows a slow increase with 
field, including a noticeable peak at higher temperature (Fig. 6). In much stronger fields 
{fioH 3> 1 T) the adiabatic response is completely suppressed as would be expected (Fig. 
6), but at an intermediate field {/hqH ^ 0.920(8) T), Xs{H) exhibits a striking peak very 
near to the (internal) field of the zero temperature phase transition. At this field the ice rule 
is locally broken and 1/4 of the spins in the sample may then be fiipped at zero energy 
cost. However, in contrast to the zero field result, Xs{T) measured near this crossover field 
(0.86 T) exhibits a simple Curie law, xs — C /T (Fig. 5), indicating a different type of 
magnetic current to that observed in the weak field limit, as anticipated in Ref.^°. We may 
regard the magnetic response in this regime as characteristic of switching magnetic dipoles, 
rather than magnetic monopoles. Note that the temperature evolution was measured at 
this point just off the peak maximum as it was found that systematic errors in fitting to 
the Cole-Cole function are minimised at this point (Supplementary Information, Section 
5). The Curie constant C may be calculated under the assumption that 1/4 of the spins 
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are thermally active and that these have a projection of 1/3 of their full classical value 
on the field direction. Thus we predict C = (7/36 = 0.1097 where C ^ 3.95 is the high 
temperature Curie constant. A fit of the experimental data to the expression xs{l^oH — 0.86 
T)= tt/(T — Tc) gave a — 0.090(5), T^; = 0.4(2) in close agreement with our prediction (Fig. 

5) . 

The striking 1/T divergence and location of the peak position in the H — T plane (Fig. 

6) suggests that the adiabatic susceptibility is dominated by the classical critical point for 
monopole condensation (Fig. 1), for which the isothermal susceptibility Xt is predicted to 
diverge as 1/|T — T^l^, with 7 = 1 and Tc here. Usually the ratio Xt/xs, analogous to 
the Landau-Placzek ratio in a fluid, should diverge towards the critical point. However, our 
data for xt and xs (Fig. 5) illustrate a fairly typical paramagnetic response, as discussed 
above, with xs ^ Xt- Hence both xs and Xt diverge, but the latter is always larger, as 
required by thermodynamics (Eqn. 5). 

According to Ref.^^, the field dependence of the susceptibility should be mean field hke, 
with logarithmic corrections: hence we would expect x ^ \H — Hc\^^^^^~^ with 5 = 3. 
However we observe an exponent of 2 rather than 2/3 (Supplementary Information, Fig. 
S7). Thus, defining reduced variables t — Tc (with Tc = 0.36 here), and h ^ H — Hc^ 
we find to a good approximation (suppressing dimensional constants): 

This implies h/^/i scahng (Fig. 6, inset), which is formally characteristic of a zero dimen- 
sional phase transition (Supplementary Information, Section 6). An alternative interpre- 
tation of the field dependence is in terms of a classical single spin fiip process, associated 
with the 'free' moments in the eventual ordered structure, which would be characterised by 
a response of the type xs ^ + ^^) and hence h/t scaling, but our data appears to 
distinctly rule against this possibility (Fig. 6, inset). Thus the behaviour of Xs{H) seems 
inconsistent with both the 'monopole' and T = oo fixed points. One possibility, as discussed 
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further in the Supplementary Information (Section 6), is that the susceptibihty is dominated 
by a zero temperature quantum critical point, but again with anomalous exponents. It is 
noteworthy that there have been reported several other examples of anomalous exponents 
in the quantum critical behaviour of rare earth magnets 

5 Conclusion 

While the concept of magnetic monopoles in spin ice is supported by much experimental 
evidence the microscopic mechanism of monopole motion has yet to be identified. Our 
investigation of Dy2Ti207 has isolated the characteristics of this mechanism to which any 
future theory must conform. Our firm result is that monopoles obey the Nernst-Einstein- 
Smoluchowski equation with temperature independent diffusion constant, a strong signature 
of Brownian diffusion. It should be emphasised that this is an experimental result and not 
a theoretical input. 

It is interesting to discuss these results in the context of band theory. Just as water ice 
can be thought of as an intrinsic protonic semiconductor^^, so spin ice can be thought of 
as an intrinsic semiconductor for magnetic monopoles. These are produced by the thermal 
unbinding (or 'fractionalisation'^) of conventional magnetic excitons. They tunnel from site 
to site and have an effective mass determined by the inverse Debye length (proportional^^ 
to ^Jx/T) . The fact that the carriers support a temperature independent diffusion constant 
places an essential constraint on any theoretical development of this picture. 

We have shown that the adiabatic susceptibility gives a new perspective on the magnetic 
properties of spin ice, revealing a direct measure of the magnetic monopole concentration 
and critical behaviour in applied field. It would be useful to apply our methods to the 
low temperature regime, as zero-field measurements in that regime await an unambiguous 
interpretation in the monopole picture and the theory of Ref.^° has yet to be compre- 
hensively tested. More generally we may conclude that the adiabatic susceptibility, often 
ignored as an uninteresting by-product of ac-susceptibility analysis, may contain a wealth of 
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information about strongly correlated spin systems at low temperature. 

Our results have revealed a new property of magnetic monopoles: their partial magnetic 
polarisation by an applied field. Taken together with the remarkable prediction that a 
monopole will carry an electric dipole moment (the equivalent of its spin, if we reverse 
the roles of electricity and magnetism) a fascinating picture of the local properties of the 
monopole is starting to emerge. In general, the properties that we have discovered will have 
an important influence on any future application of magnetic monopoles in spin ice that seek 
to exploit their local quantum degrees of freedom. 

Methods 

The dynamical magnetization of a 0.0326(1) g cubic crystal of Dy2Ti207 was measured 
with the ACMS (AC-Measurement System) option of a PPMS (Physical Property Measure- 
ment System, Quantum Design). Alternating and direct current magnetic fields {Hac and 
Hdc^ respectively) were applied parallel to the cubic [111] axis of the sample. Data were col- 
lected at different temperatures between 1.9 K < T < lA K in the ac- frequency range of 10 
Hz to 10 kHz. A variable dc field of /hqH = — 10 T was apphed (at low field the absolute 
field was calibrated in dc-sweep measurement). Scans were taken at different ac fields in 
the range ji^Hac = 0.05 — 3 x 10~^ T to dispel the possibility of non-hnear response of the 
system. The results presented here were taken at /io^ac = 5 x 10~^ T. Data were corrected 
taking into account a demagnetizing factor V — 1/3 to give Xac — ^ac/ {{Hac — T^Mac)- The 
calibrated response function of the instrument was checked by measurement of a very dilute 
paramagnetic salt (Supplementary Information, Fig. S5). 

The data were fitted to the phenomenological model for the frequency dependent sus- 
ceptibility described in Ref.^^'^^. By separating Eqn. 1 it is possible to derive analytical 
expressions for the real and imaginary parts and argand diagram (Cole-Cole plot), which 
were each fitted to the experimental data at a given temperature using a single set of pa- 
rameters xs, XT, T, a. 
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The Cole-Cole formalism assumes a symmetric unimodal distribution of logarithmic re- 
laxation times Inr' with mean ln(T). If the variance a'^^^, <^ 1 then it closely approximates 
al where z/ = 1/r' is the relaxation rate. It can be shown that^^: 



In the text we label crf^^, simply as crf^^ for ease of reading, although strictly r is a fixed 
parameter at a given temperature. It should be noted that there is no general way to derive 
the true mean relaxation time {r') from ac-susceptibility data: here we approximate it to 
the Cole-Cole parameter r. 

The dimensionless monopole density x{T) was estimated by fitting experimental specific 
heat data to Debye-Hiickel theory The specific heat was represented as the temperature 
derivative of the energy per diamond lattice site: 



where /i is the monopole chemical potential, jJiDniT) is the Debye-Hiickel correction, cal- 
culated self consistently with the dimensionless monopole density x(T), and udcm{T) is a 
correction for double charge monopoles. The experimental specific heat data, taken between 
0.4 K < T < were fitted by adjusting /i, with the best fit value ji/k — —4.33 K. 

The theory is not exact in the temperature range of interest and Figs. 4 and 5 report an 
approximate envelope of systematic error, found by extrapolating the theory between low 
and high temperature according to different schemes (Supplementary Information, Section 



To derive Eqn. 4, assume r ^ oc i7, take logarithms and Taylor expand the right hand 
side of Eqn. 3 to find Inr = /ii + h2\[x\ then if h\ and /i2 are uncorrelated, Eqn. 4 follows. 




(8) 



(9) 



7). 
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Figure 1: Right: Fragment of the crystal structure of spin ice, illustrating the 
hopping of emergent magnetic monopoles. Fragment of spin ice's cubic pyrochlore lat- 
tice, which consists of corner-hnked tetrahedra, showing spin configurations (arrows). Top 
panel illustrates crystallographic axes, the applied field direction and how internal fields 
may be transverse to the local spin direction. Blue (red) circles represent negative (posi- 
tive) monopoles. Bottom panel illustrates how a monopole hop can be associated with a 
spin fiipped by a transverse field or a tunnelling event through a potential barrier. Left: 
Temperature (T) versus field {R || [111]) phase diagram of spin ice Dy2Ti207. The 
full line is a line of first order phase transitions, terminating in a classical critical point, that 
has been interpreted as a monopole condensation^. Monopoles are deconfined in zero field 
but become confined in an applied field. The right hand diagrams show how the monopoles 
reform fiippable spins or dipole pairs near the critical field. Dotted lines are guides to the 
eye. Experimental points with error bars show the applied field of the maximum in the 
adiabatic susceptibility measured here. 
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Figure 2: Experimental observation of a finite adiabatic susceptibifity (xs) and its 
relation to the isothermal susceptibility (xt)- These quantities are estimated as the 
real part of the frequency dependent susceptibihty x(^) the hmits uj ^ oo and ^ 0, 
respectively as indicated in the main plot. Experimental data are at T = 1.95 X at applied 
static magnetic field /xo|i/| = (black circles), 0.86 (grey circles) and 10 T (light grey circles). 
At 10 T both susceptibilities are nearly zero. At 0.86 T, the two are of similar magnitude. 
At zero applied field Xs ^ Xtj but still finite. The respective lines are the fit to a Cole-Cole 
model (see text). The inset shows the /hqH = data on log-log scales. The clear deviation 
from a linear curve at large uj confirms the presence of a finite offset, xs- Here the blue line 
is the fit to the Cole-Cole function using finite xs and the red line is the same fit with xs 
constrained to be zero. 



21 
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Figure 3: Ac-magnetricity and monopole diffusion in spin ice, Dy2Ti207. Left: 
a representative fit to the ac-magnetization data at Hdc = 0: (middle) real and (lower) 
imaginary susceptibilities, and (upper) Cole-Cole plot (an argand diagram). The lines (red, 
T = 4.5 K; blue, T = 8 K; green, T = 14 K) fit all data with the same set of four parameters at 
a given temperature: Xt, Xs^ <^ and r, as defined in the text. Right (main figure): monopole 
mobility measured at applied fields /io|^| = (full black circles), 3 (full dark grey), 10 
(full light grey), 18.5 (open black) and 38.5 (open grey) mT. The red line is u = A/T (red 
line) where A = 2.27(2) x 10~^ m K s~-^ T~-^. This is characteristic of Brownian diffusion 
of monopoles with a temperature-independent diffusion constant (inset). The blue line is 
u = Be~^/^ (blue line) where B = 39(1) m s"-*^ T"-*^, C = 250(1) X, characteristic of 
a previously identified Orbach-like spin fiip process arising from the coupling to excited 
crystal field states. This process is extinct below 10 X, giving way to monopole diffusion as 
the cause of magnetic relaxation in spin ice. 
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Figure 4: Monopole signatures. The experimental variance in logarithmic relaxation 
time a^^^ (circles) compared with the predicted form for monopolar-field assisted tunnelling. 
Red and green indicate, respectively, applied fields oi fioH = and the set of finite fields listed 
in the caption of Fig. 3. For each curve there are two fitted parameters and a|, which 
describe the mean square static field and mean square monopole field, respectively. The line 
is the function + xa^ where x{T) is the monopole density. The shading indicates the 
maximum possible systematic error in the monopole density. We find ai^2 = 0.84(1), 3.40(5) 
respectively in zero field and 1.15(1), 3.40(5) in finite field. The deviation at T > 10 X is 
related to a change in relaxation mechanism (see text). 
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Figure 5: Monopole and dipole signatures revealed by comparing the measured 
adiabatic {xs) ^^id isothermal (xt) susceptibilities. Upper figure: the susceptibilities 
are shown as a function of apphed dc-field at zero field (black circles) and near the crossover 
field (grey circles) shown as a red dotted line in Fig. 1. The full red line is the experimentally 
measured monopole density fitted to Xs'(O) by the adjustment of a scale factor (see text). 
The blue line is a Curie- Weiss law fitted to Xt(0) (Supplementary Information, Fig. S4). 
The light blue line is the theoretical prediction xs{Hc) = C /{T — Tc) (see text). Lower 
figure: Measured xs as a function of weak applied field < /xgi/ < 0.05 T (circles, same 
colour code is maintained as in Fig. 3). The full red line is the experimentally measured 
monopole density fitted to Xs{^) by the adjustment of a scale factor (see text), also showing 
the maximum systematic uncertainty in the measured monopole density (red shading). 
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Figure 6: Evidence of unconventional critical behaviour. Main figure: adiabatic 
susceptibility versus applied magnetic field, showing an unusual Lorentzian field dependence, 
suggestive of anomalous critical behaviour. The figure illustrates zero response at strong field 
and a peak at the field where the ice rule breaks. The lines are fits to a Lorentzian function 
at the following temperatures: blue, T = 1.95 K; light blue, T = 2.5 K; purple, T = 3.0 
K; lilac, T = 4.0 K; light green, T = 6.0 K; green, T = 8.0 K. Inset: scaling of Lorentzian 
width (A): the red line is the function A oc Experimental data appears to rule out the 
possibility of a diflFerent scaling - e.g. linear scale h/t. 
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Brownian Motion and Quantum Dynamics of Magnetic 
Monopoles in Spin Ice: Supplementary Information 



L. Bovo^, J. A. Bloxsom^, D. Prabhakaran^, G. Aeppli^, S. T. Bramwell^. 
1. London Centre for Nanotechnology and Department of Physics and Astronomy^ 
University College London^ 17-19 Cordon Street, London, WCIH OAH, U.K. 
2. Department of Physics, Clarendon Laboratory, University of Oxford, Park Road, Oxford, 

0X1 3PU, U.K. 

1. Model Used to Fit the data 

In the main text we note that the relaxation time r does not have the same physical in- 
terpretation as that arising in conventional models of paramagnetic spin relaxation. However 
this fact may be neglected for fitting purposes, and the data analysis reduces to a familiar 
problem of fitting x{^) to conventional phenomenological forms. 

The Debye model, with a single exponential relaxation of the susceptibility fails in the 
description of experimental data. Various empirical equations have been formulated to give 
the variation of Im[x(ct;)] = ^' with Re[x(cc;)] = ^ as the frequency is varied, taking into 
account the presence of different distributions of relaxation times. In particular, as discussed 
in the main text, we tested two of the most commonly used models: the Cole-Cole (CC)^ 
and Cole-Davidson (CD)^ formahsms. Fig. SI shows an illustrative comparison between the 
two representations. The CC formalism works rather well in all range of frequency for a wide 
interval of temperature, 3.5 <T <1A K. At low temperatures, in the low frequency regime, 
the experimental data start to deviate form the model, showing an increasing asymmetry. 

Nevertheless, it is possible to fit all data with the same set of four parameters for each 
temperature (red lines in Fig. SI show two examples for data taken at T = 1.95 K and 6 K, 
respectively). On the contrary, when the CD-model is applied, the fit to the experimental 
data is poorer, as shown in Fig. SI. In particular, it is not possible to find a unique set 




Figure SI: Comparison between Cole-Cole (CC) and Cole-Davidson (CD) model. 

Top) Data taken at zero dc-field and T = 1.95 K: (upper small) real and (lower small) 
imaginary susceptibility and (main) argand diagram. Bottom) Data taken at zero dc-field 
and T = 6 K: (upper small) real and (lower small) imaginary susceptibility and (main) argand 
diagram. CC model: The red line fits all the data with the same set of four parameters: 
XtjXs^^ and r, as defined in the main text. The deviation that happen at low temperature 
and low frequency is discussed in the text. CD model: The blue and turquoise lines represent 
the two best fit obtained analysing experimental data. It must be noted that, unlike the 
CC model, there is not a unique set of parameters that can fit simultaneously the three 
diagrams. 
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of parameters to simultaneously fit the real and imaginary frequency dependence of the 
susceptibility as well as the argand diagram {xt ^ fiXr))- 1^ clear then that the CC 
formalism seems to give a more representative picture of the dynamics in Dy2Ti207. This 
seems physically reasonable as the CC model assumes a distribution of logarithmic relaxation 
times that is cut off exponentially at high frequency, while the DC model assumes a power 
law to high frequency. The latter is not consistent with what is known about the magnetic 
dynamics in spin ice, where the monopole hop rate would presumably lead to a cut-off at 
high frequency. Hence the superior performance of the CC model is plausible. Undoubtedly, 
the determination of the proper distribution of relaxation times for a spin-ice system could 
in principle lead to a better interpretation of experimental data, especially of the asymmetry 
showed at low temperatures in the low frequency regime. 

2. Relaxation Time 

The magnetic relaxation time as a function of temperature of Dy2Ti207 is displayed in 
Fig. S2. The general dynamic behaviour, in good agreement with previous experimental 
data^, is consistent with the presence of two different regimes. At enough high temperature, 
above 10 the time scale drops dramatically due to a thermally activated process (Orbach- 
like mechanism) with an energy barrier of few hundreds kelvin compatible with higher energy 
crystal field levels being populated. This barrier, Ea ~ 250 K, is insurmountable below 10 

and the system enters a quasi-plateau region. A further lowering of the temperature, 
below 2 K - outside the temperature range investigated in this work - would have caused a 
sharp upturn of the relaxation time^. Following the work done by Jaubert and Holdsworth^, 
it is possible to interpret the quasi-plateau region with the presence of a thermal assisted 
quantum tunnelling through the crystal field barrier: below 10 K the spins are Ising like and 
the system can be represented by stochastic single spin dynamics, or in 'monopole' language 
by the creation and diffusion of the topological defects. 
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Figure S2: Temperature dependence of relaxation time. Experimental data ex- 
trapolated from the CC-fitting (black circles). The red line is a fit to the free diffusion of 
topological defects in the nearest neighbour approximation^. Inset: same data shown on 
a different scale. The blue line represents the fit to an Arrhenius law, with an activation 
energy of about 250 compatible with a classical Orbach-like mechanism (main text). 

In a first approximation^, it is possible to consider an Arrhenius law of the type: 

r = To exp(2 Je///fcr), (1) 

where tq is the microscopic tunnelling time and ^J^ff is the energy cost of a single free 
topological defect in the nearest neighbour approximation and is half of that of a single spin 
fiip. The red line in Fig. S2 represents the best fit to Eqn. 1 with tq = 8.6(2) x 10""^ s and 
2Jeff/k = 2.8(1) K, in close agreement with the results of Jaubert and Holds wort h^ and 
in line with the analysis and description that we have proposed in the main paper. This 
expression has been proven to work rather well in the semi-plateau regime, but it starts to 
fail below 2 K underestimating the relaxation time at very low temperature. Nevertheless, 
they showed how it is possible to interpret the magnetic relaxation of Dy2Ti207, in the whole 
temperature range in terms of the diffusive motion of monopoles in the canonical ensemble, 
constrained by a network of 'Dirac strings' filling the quasi-particle vacuum. 

A similar way to rephrase the problem is to consider that at any given temperature there 
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is only a well defined fraction x of fiippable spins (or density of monopoles, reverting to the 
topological defects representation) that can directly contribute to the magnetic relaxation 
process. Spins that are not associated with monopoles can flip at much hinger cost and 
cannot contribute to the dynamics. In this work, indeed we have shown that once taken into 
account the thermal evolution of the density of monopoles and the temperature factor char- 
acteristic of a diffusive (Brownian) motion (main text), the resulting characteristic hopping 
rate turns out to be temperature independent, symptomatic of a spin relaxation that occurs 
by quantum tunnelling. 

3. Ambiguities in the Interpretation of the Dispersion of Relaxation Times 

It should be noted that Eqn. 4 of the main text is only strictly valid for /ii, /i2 ^ 1 which 
is marginal for the data considered here. Unfortunately the inversion of a logarithmic time 
distribution to a frequency distribution is an ill defined mathematical problem^ that can 
only be accomplished in a satisfactory way for narrow distributions. This means that our 
simple model is not ruled out, but nor is it unambiguously implied by our data. In contrast 
an unambiguous deduction from our data is that the variance of effective energy barriers is: 

al = {kTf{al + xal), (2) 

where cri,(j2 now become uninterpreted parameters. This quantity is plotted in Fig. S3. 

While Eqn. 2 is weaker than Eqn. 4 (main text) at the level of physical interpretation, 
the presence of factors of kT in its right hand side clearly imply a quantum relaxation and 
its factor of x implies cooperative dynamics involving magnetic monopoles. The direct field 
interaction that we have suggested is one plausible explanation of this. 

A minimal conclusion that may be drawn from our analysis is as follows. Assuming 
the monopole model, we can directly measure the diffusion constant D{T). However, this 
argument cannot be reversed, so the behaviour of the measured relaxation time using the 
Cole-Cole model does not in itself infer the role of magnetic monopoles. This ambiguity may 
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Figure S3: Variance of the effective energy activation barrier versus temperature. 

Temperature dependence of the experimental variance in energy scale as a function of weak 
applied field < fioH < 0.05 T (circles, same colour code is maintained as in Fig. 3 of the 
main text). Red and green indicate, respectively, applied fields of fioH = and the set of 
finite fields. The line is Eqn. 2 with ai^2 = 0.84(1), 3.40(5) respectively in zero field and 
1.15(1), 3.40(5) in finite field (Fig. 4 main text). The deviation at T > 10 K is related to a 
change in relaxation mechanism (main text). 

be traced to the fact that the relaxation due to monopoles takes the form of an effective spin 
relaxation (see above). In contrast, the fact that the dispersion of the logarithmic relax- 
ation time correlates with the measured monopole density does directly implicate magnetic 
monopoles in the relaxation process. 

4. Isothermal Susceptibility 

For the sake of completeness. Fig. S4 reports the evolution of isothermal susceptibility as 
a function of temperature in the limit of weak applied field {fioH < 0.05 T). For comparison, 
the bulk magnetometry SQUID data (blue circles) of the same sample are also displayed in 
Fig. S4. The best fit (red line) to an apparent Curie- Weiss law Xt = Xc/{T — 9) gives rise 
to a Curie Constant xc = 4.25(5) K and 9 = 0.8(1) K. Besides, experimental data deviate 
from this fitting at low temperature: the susceptibility does not approach a constant value 
as T approaches zero. Similar behaviour was reported for another spin ice compound^'''. As 
already mentioned in the main text, our results are in good agreement with the prediction^ 
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Figure S4: Isothermal susceptibility. Results obtained fitting experimental data in the 
limit of weak applied field (0 < /hqH < 0.05 T); the same colour code is maintained, as in 
the main paper. The blue circles represent the bulk SQUID magnetometry measurement. 
The red line is the fit to a Curie- Weiss behaviour as reported in the text. 

that, although the crossover of the Curie constant for a spin ice system is expected to be 
incomplete in the temperature range studied, the experimental value is higher than the 
expected xc = 3.95 K for Dy2Ti207. 

5. Adiabatic Susceptibility 

To further assess the consistency of our fits, we tested the reliability of the extrapolation 
of xX^) to the limit of ^ oc, to give xs- First to exclude the possibility of the presence 
of an apparent background arising from the instrumental response, we determined the latter 
by measuring a single crystal of paramagnetic LiYF4 doped with 0.3 at. % Ho. As shown in 
the right panels of Fig. S5, the signal is at least two orders of magnitude lower than in the 
case of the spin-ice crystal in the entire frequency and temperature range. Furthermore, as 
shown in the left panel of Fig. S5, the temperature dependence of the experimental value 
of x^ for Dy2Ti207 measured at the highest applied frequency {uj = 10^ s~^) replicates the 
temperatures dependence of the estimated xs{T)^ showing that xX^) tends analytically to 
the high frequency limit. 

As described in the main text, we also measured the ac-susceptibility as a function of 
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Figure S5: Left. Experimental adiabatic susceptibility: (grey circles) extrapolated 
from CC-fitting and (black circles) the measured value of xX^O^ ^~^)' The red line represent 
the expected monopole density adjusted by a scale factor; the shaded area indicates the 
maximum systematic error in the estimated monopole density (main text and Fig.s 4 and 
5). The blue line, and respective shadow, is obtained simply by shifting the red line of an 
offset of about 0.0105(1). Right. Instrumental response: Real (upper) and imaginary 
(lower) component of the susceptibility measured at different temperatures for a single crystal 
of LiYF4 : 0.3 at. % Ho; the signal is at least two orders of magnitude lower than in the case 
of the spin-ice crystal in the entire temperature range. 




X" X' 

Figure S6: Cole-Cole plot and respective CC-fit of data taken at T = 1.95 K at 
different applied field. Right: representative data taken at, from dark to light green, /lcqH 
of 0.1, 0.2, 0.3 and 0.4 T respectively. Left: representative data taken at, from dark blue to 
turquoise, /hqH of 0.5, 0.7, 0.8 and 0.9 T, respectively. 
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Figure S7: Adiabatic susceptibility versus applied field, Log-Log plot. Experi- 
mental xs extrapolated from data at T = 1.95 K. The blue line shows that the observed 
exponent of a = 2 for xs ^ \H — Hc\^~^'^ • 

an applied static magnetic field along [111]; Fig. S6 shows argand diagrams obtained at 
T = 1.95 K. First of all, it is worth notice that the deviation in the low frequency regime 
of the fitting from the experimental behaviour decreases with the increase of the applied 
field. Because of the low signal measured at high field, the complete analysis was performed 
only on data taken at jUoH < 1 T; at higher field only the frequency dependence of x\^) 
was fitted. For this reason, we preferred to measure the evolution of Xs{T) as a function 
of temperature at fioH = 0.86 T instead of the actual crossover field where xs showed its 
maximum value (main text and Fig.s 5 and 6). 

Looking at the field dependence of the adiabatic susceptibility, according to mean field 
theory^, one would expect xs ^ \H — i^d^"^/^^. Fig. 7 shows as an example data taken at 
T = 1.95K. The Log-Log plot clearly indicates an exponent of 2 rather than 2/3. 

6. Critical Behaviour 

If we interpret our field dependent measurements in terms of classical critical exponents, 
then we are led to the unusual conclusion that the critical point involved is zero dimensional. 
Thus, if we assume a zero dimensional critical point {d = 0) then the scaling laws predict 
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the set of exponents 6 — —1^ j — 1^ /3 — —1/2 and ly and rj are of course undefined. Here 
we show how these numbers summarise the experimental observations. 

In critical point theory we can describe the field and temperature dependence of the 
order parameter m by 'polar' co-ordinates r, 9^ where r is related to the distance from the 
critical point in the {h, t} plane and 9 is related to the 'angle' with respect to the t axis. 
Here t ^ T — Tc, h ^ H — He, m ^ M — Mc- The polar equations in their most general 
form: 

h^T^'Ue), (3) 
t^ftir,e), (4) 

m ~ r^Uie), (5) 

where the three functions ft,m,h on the right hand side are to be determined. 
In mean field theory the transformation can be written expUcitly: 



h ~ r^'eii - e') (6) 

i~r(l-2^2), (7) 
m - r^e, (8) 



with the set of exponents 7 = l,/3 = l/2,(5 = 3. 
If we write the equations, 



h-^r^^i&ud, (9) 
t - r, (10) 



m ^ r 



^0, (11) 



with the exponents 7 = !,/? = — 1/2, 5 = — 1, we recover our result x = l/{t + h'^). The 
negative critical exponents look unusual but obey thermodynamics and scaling theory, as is 
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evident from the preceding equations. For example, the Griffiths scahng relation is obeyed: 
7 = — 1) and m is always an increasing (arctangent) function of h. Applying the usual 
scaling relations we find: a — 2^ d — and 77, u are undefined as one would expect for a zero 
dimensional critical point. As the correlation length and its exponent u are undefined, the 
dynamical exponent z is also undefined. 

Alternatively we could consider a single classical spin in a magnetic field h. It is easy to 
show that to small h the susceptibility is approximated by x ^ T/ [h? -\-T'^). In this case the 
polar equations can be satisfied as above, but with /3 = 5 = 0; however, the Griffiths scaling 
relation is not obeyed, suggesting that the classical single spin problem cannot be formally 
represented as a critical system in the zero temperature limit. 

A final possibility is to interpret the exponents in terms of a quantum critical point. Thus, 
although the phase transition at zero temperature is first order, at a significant 'distance' 
away the appropriate renormalisation fiow might be dominated by a T = critical point 
of a similar sort to that which occurs the 'transverse field Ising model', but which is in 
reality removed by the Coulomb interaction, a strongly relevant variable. A quantum critical 
point typically exhibits the classical critical exponents appropriate to one higher dimension^. 
For a three dimensional system this means that quantum critical points are at the upper 
critical dimensionality of 4, and so exhibit mean field exponents. As shown above, the 
observed 7 exponent is mean field like, but when combined with the apparent 5 exponent of 
5 — —1 suggests a different universality class. In the theory of quantum critical points^, this 
change in dimensionality could be accommodated by a negative dynamical [z) exponent, but 
this would be physically hard to justify. Thus the field dependent behaviour is distinctly 
anomalous, however it is viewed. 

7. Specific Heat 

The specific heat of single crystal Dy2Ti207 was measured with a Quantum Design Physi- 
cal Properties Measurement (PPMS) System. The Debye-Hiickel theory of Ref.^^, along with 
a correction for double charge monopoles^^ was used to fit the specific heat data over the 
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entire range of measured temperature, and so to extract the monopole charge density as a 
function of temperature. The monopole chemical potential was used as a fitting parameter 
and was refined to be ^ 4.33 K, which is in good agreement with theory The Debye-Hiickel 
theory is accurate only for small x/T. This means that it is accurate at low temperature, 
and also a reasonable approximation at high temperature, but breaks down at intermediate 
temperature, near to the peak in the specific heat at ^ 1 iC. Although this leads to a 
significant systematic error on the specific heat, the error that propagates through to x{T) 
is relatively small. Estimates of this systematic error were made by altering parameters in 
the partition function within reasonable bounds, to give the envelope of curves exhibited in 
the paper. The estimated x{T) was found to be in close agreement with that derived from 
numerical simulations^^ and is clearly sufficiently accurate for our purpose. Full details of 
the specific heat analysis will be published separately. 
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